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ABSTRACT 



Context. The evolution of a young protocluster depends on the relative spatial distribution and dynamics of both stars 
and gas. 

Aims. We study the distribution and properties of the gas and stars surrounding the luminous (10^ Lq) protocluster 
IRAS 18511+0146. 

Methods. IRAS 18511+0146 and the cluster associated with it has been investigated using the sub-millimetre (JCMT- 
SCUBA), infrared (Spitzer-MIPSGAL, Spitzer-GLIMPSE, Palomar) and radio (VLA) continuum data. Cluster simula- 
tions have been carried out in order to understand the properties of clusters as well as to compare with the observations. 
Results. The central most obscured part of the protocluster coincident with the compact sub-millimetre source found 
with SCUBA is responsible for at least 2/3 of the total luminosity. A number of cluster members have been identified 
which are bright in mid infrared and show rising (near to mid infrared) spectral energy distributions suggesting that 
these are very young stellar sources. In the mid infrared 8.0 fim image, a number of filamentary structures and clumps 
are detected in the vicinity of IRAS 18511+0146. 

Conclusions. Based on the luminosity and cluster size as well as on the evolutionary stages of the cluster members, IRAS 
18511+0146 is likely to be protocluster with the most massive object being a precursor to a Herbig type star. 



Key words. Stars 
18511+0146 



formation - Stars: pre-main sequence - Infrared: ISM - Submillimeter - Stars : individual : IRAS 



1. Introduction 

Massive stars form accompanied by swarms of lower mass 
objects and the relationship between the two is by no means 
clear. It is however known that the development of a young 
cluster depends sensitively on events which occur in the 
phase when it is still surround ed by the re mnants of the 
"core" from which it formed (jLada fc Lad a 2003). It is 
therefore of interest to obtain sensitive observations of the 
young protostars and the surrounding gas in the phase 
when dust obscuration causes the cluster to be invisible at 
optical wavelengths. The objects of interest are thus best 
observed at infrared and radio wavelengths. 

This article reports such a study of the cluster sur- 
rounding what appears to be a young intermediate (2-8 
Mq) protostar (IRAS 18511+0146, also known as Mol 75 
and RAFGL 5542) from the survey of lMolinari et al.l (|l996l 
11998. ) using both data from the Spitzer satellite and ground 
based observations. We proceed on the hypothesis that the 
cluster surrounding IRAS 18511+0146 is essentially a fore- 
runner of the small clusters examined by Testi et al. (1998). 
Thus, the bolometric luminosity (roughly 10"* Lq) and size 
(0.5 parsec) are similar but, as discussed later, the visual 
extinction is large (of order 50 visual magnitudes) and vari- 
able. 
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The gas distribution is of great importance both because 
of its effect upon the extinction to individual protostars 
and because the gas plays an important role dynamically. 
The outflows and ionization caused by the young protostars 
are thought to be the primary agent causing the cluster to 
disperse. We therefore use JCMT-SCUBA observations to 
assess the mass and distribution of the gas associated with 
IRAS 18511+0146. We then use this information in simu- 
lations which we have carried out with the aim of testing 
our hypothesis that the young protocluster associated with 
IRAS 18511+0146 is an obs cured ve rsion of the intermedi- 
ate mass clusters studied by iTesti et al.i (jl998i ) . 

Previous studies o f IRAS 18511 were carried out by 
IWatt fc Mundvl (|1999D using the OVRO interferometer and 
the VLA to map both molecular line and continuum emis- 
sion. They detected what appeared to be a small ionized 
region ionized by a Bl sta r which coin c ided with a com- 
pact clump seen in C^^O. iKurtz et al.l (|2004j) detected a 
methanol maser which was however offset by roughly 0.4 pc 
(^^19") relative to the ionized gas fo r a distance of 3.9 kpc 
to IRAS 18511 (jMolinari et al.lll996f ). lMiii et al.l (|2002D de- 
tected the 3 /zm water ice feature towards the central object 
suggesting a visual extinction o f 22 magnitudes along this 
line of sight. [Zhang et ahl (|2005f ) have reported strong indi- 
cations of an outflow in IRAS 18511 using the ^■^CO line. 
According to them, the mass- loss rate is ^3. 7xl0~^ Mq 
yr~^ and the dynamical timescale is 5x10"^ vrs. lBrand et al.l 
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(|2001h imaged IRAS 18511 using the IRAM 30-m and 
KOSMA telescopes in a number of molecular lines (^^CO, 
CS, HCO+) and find blue shifted and red shifted compo- 
nents distributed over various clumps. 

In this study, we explore in detail the region sur- 
rounding the IRAS source 18511-1-0146 using infrared, sub- 
millimetre, and radio continuum data. In Sect. 2, we sum- 
marize the observations. In Sect. 3, we present the observa- 
tional results and in Sect. 4, a model of the relative distri- 
butions of gas and stars is examined. A discussion of IRAS 
18511 region is presented in Sect. 5, and a short summary 
of our conclusions is given in Sect. 6. 

2. Observations, Available data and Data reduction 

2.1. JCMT-SCUBA 

The sub-millimetre observations of IRAS 18511 at 450 and 
850 /im, using the Submillimctrc Common User Bolometer 
Array (SCUBA) of the James Clerk Maxwell TelescopeQ 
were carried out on 28 May 1999. The data were processed 
using their standard pipeline SCUBA User Reduction 
Facility (SURF). The planet Uranus was used for calibra- 
tion. Submillimetre maps were generated at 450 and 850 /im 
and the fluxes extracted. The beam size is 10" at 450 /im 
and 15'.'5 at 850 fim. The sensitivities are 0.3 Jy/beam 
and 0.04 Jy/beam at 450 and 850 ^m, respectively. The 
flux densities extracted were used to construct the Spectral 
Energy Distribution (SED) of IRAS 18511. 

2.2. Spitzer data 

The Spitzer Space TelescopeEl (| Werner et all 120041 ) was 
launched in space in August 2003 and consists of a 
0.85-meter telescope with three cryogenically cooled in- 
struments: InfraRed Array Camera (IRAC), Multiband 
Imaging Photometer for Spitzer (MIPS) and InfraRed 
Spectrograph. We have used data from the Legacy projects, 
GLIMPSE (PI: E. Churchwell) and MIPSGAL (PI: S. 
Carey) in this paper. The GLIMPSE and MIPSGAL im- 
ages have been obtained using the software 'Leopard'. 

2.2.1. MIPSGAL 

The Multiband Imaging Photometer for Spitzer (MIPS) 
provides the Spitzer Space Telescope with capabilities for 
imaging and photometry in broad spectral bands centred 
nominally at 24, 70, and 160 /im, an d for low-resolutio n 
spectroscopy between 55 and 95 /im (Rieke et al.' '200j). 
The Multiband Imaging Photometer for Spitzer Galactic 
Planey survey (MIPSGAL) surveyed the sky in identical 
regions as covered by GLIMPSE (nex t subsection) at 2 4 
and 70 ^J.m using the MIPS instrument (jCarev et al.il2005D . 
The instrument achieves diffraction-limited resolution of 6" 
and 18" at 24 and 70 /im, respectively. The pixel size is 

^ This paper makes use of data from the James Clerk Maxwell 
Telescope (JCMT) Archive. The JCMT is operated by the Joint 
Astronomy Centre on behalf of the UK Particle Physics and 
Astronomy Research Council, the National Research Council of 
Canada and the Netherlands Organisation for Pure Research. 

^ This work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, under NASA contract 1407. 



2'.'55 at 24 /im and 9"98 at 70 /im. The Spitzer Science 
Center provides final mosaics called Post-Basic Calibrated 
Data (PBCD) products. These are the maps of multiple 
calibrated image frames or BCDs (Basic Calibrated Data). 
For IRAS 18511, we found that these products satisfied our 
scientific goals. We have carried out extraction and photom- 
etry of the sources in the region around IRAS 18511 using 
the MIPSGAL PBCD images at 24 /im and 70 /im (details 
given in Appendix A). The details of the extracted sources 
and their fluxes are presented in Section 3. 

2.2.2. GLIMPSE 

In the GLIMPSE (G alactic Legacy Infrar ed Midplane 
Survey Extraordinaire: [Benjamin et ahl [20031 ) project, the 
Spitzer Space Telescope surveyed approximately 220 square 
degrees of the Galactic plane covering a latitude range 
of \b\ < 1° and a longitude range of 10° < / < 65°, 
-65° <l < -10°. This survey is carried out in the 4 IRAC 
bands. IRAC is a four-channel camera that provides simul- 
taneous 5.'2x5f2 images at 3.6, 4.5, 5.8, and 8.0 /im with 
a pixel size of l'.'2xT.'2 ([Fazio et al.ll2004f ). The resolutions 
achieved by IRAC are 2"4, 2"4, 2"8 and 3"0 in the 3.6, 4.5, 
5.8 and 8.0 /im bands, respectively. 

The sources in the IRAS 18511 region have been ex- 
tracted from the GLIMPSE More Complete Archive. The 
GLIMPSE archive contains point sources with peak signal- 
to-noise ratio greater than 5 in at least one band. The ex- 
tracted sources have been used in constructing the spectral 
energy distributions. However, the bright/saturated sources 
are not extracted in the GLIMPSE-catalog. In order to 
get a lower limit on the fluxes of these sources, we have 
carried out aperture photometry of these bright /saturated 
sources on GLIMPSE-IRAC PBCD images. For these 
bright sources, we have taken an aperture radius of 10 pix- 
els and a sky annulus of 15 pixels in order to estimate the 
fluxes. The approximate centre was determined from the ra- 
dial profiles using the task 'imexamine' in Image Reduction 
and Analysis Facility (IRAF). For such a combination of 
aperture radius and sky annulus, no aperture correction 
needs to be applied. 

The PBCD images have been used to study the spatial 
distribution of sources as well as near and mid infrared 
emission from this region. 

2.3. Palomar Data 

Deep near infrared observations of a ~3'x3'field surround- 
ing IRAS 18511 were obtained on 26 July 1999 with the 
Palomar Observatory 60-inch telescope equipped with the 
near infrared IRC-NICM0S3 camera. The field was ob- 
served in three broad bands: J, H, and Kg. A standard 
dithering technique was used to efficiently remove the sky 
emission and correct for hot or dead pixels in the 256 x 256 
NICM0S3 detector. The data were processed using recipes 
in IRAF to produce the final images. Photometric calibra- 
tion was ac hieved by ob s erving a set of standard stars from 
the lists of [Hunt et all ([1998| ) and iPersson et all ([19981 ). 
Astrometric calibration was performed tying the observed 
positions of bright and isolated sources in the Palomar field 
to the corresponding entries in the Two Micron All-Sky 
Survey (2MASS) database. We estimate this procedure to 
be accurate within 0"5. 



S. Vig, et al.: IRAS 18511+0146: a proto Herbig Ae/Be cluster? 



3 



The limiting magnitudes of our observations are found 
to be 18.1, 18.0, and 16.9 in the J, H, and bands, re- 
spectively. 

2.4. VLA 

We have also analysed the data for IRAS 18511 at 8.5 and 
15 GHz taken from the NRAO Data Archive (project ID: 
AD406) obtained using the Very Large Array (VLA) in 
the D configur ation. The ob s ervati on was carried out in 
Dec 1997 by W att fc Mu^ (|1999D . The total on-source 
integration time in each band is 30 min. The flux cali- 
brators were 1331+305 and 0137+331. 1832-105 was used 
as the phase calibrator. The NRAO Astronomical Image 
Processing System (AIPS) was used for reduction of the 
data. The beam sizes are 8'.'5x7'.'2 and 5'.'2x4'.'4 at 8.5 and 
15 GHz, respectively. The rms noise in the maps are 0.05 
mJy/beam and 0.13 mJy/beam at 8.5 and 15 GHz, respec- 
tively. 

3. Results 

3.1. Submillimetre emission from cold dust 

The JCMT-SCUBA maps at 450 and 850 /im trace the 
emission from cold dust in the region of IRAS 18511. 
The JCMT maps are shown in Fig. [T] The integrated 
flux densities upto 10% contour levels are 139 Jy and 14 
Jy at 450 and 850 /im, respectively. The peak flux den- 
sities are 1.33 Jy/beam and 2.28 Jy/beam at 450 and 
850 /im, respectively. The half-power sizes are 0.4 x 0.6 
pc^ at 450 fj,m. and 0.5 x 0.7 pc^ at 850 /im. The emis- 
sion maps show the presence of a core (q;j2ooo, ^J20oo ~ 
18'' 53"* 38!00, +01°50'29'.'l) as wefl as extended emission, 
particularly to the south. The cold dust emission shows 
an extension towards the south-east (q;j20oo, <5j20oo ^ 18'* 
53"^ 38!13, +01°50'21'.'6) which veers towards west further 
south. Emission towards the north-east of the core is also 
discerned. 

The emission at 450 fim as well as at 850 /im can be used 
to estimate the dust temperature of the core as follows. We 
know that the flux density, F^, for optically thin emission 
can be written as 

= nB,{Td)T, (1) 

where VI is the beam solid angle, Bi, is the Planck function, 
Ti, is the optical depth at frequency v and Td is the dust 
temperature. Assuming, r^^ oc , the ratio of flux densi- 
ties at any two wavelengths is a function of Td and /3. We 
have assumed (3 = 2 and used the ratio of fluxes (after con- 
volving the 450 /im map to the beam resolution of 15'.'5 
corresponding to the 850 /im beam) in order to obtain a 
dust temperature of 37 ± 5 K at the peak of the emission. 
The total m ass of the cloud ha s been estimated using the 
formalism of lHildebrandl ()1983D . Assuming an average dust 
temperature of 20 — 30 K for the entire cloud and using the 
value of dust emissivity at 850 /i m (Kssnum — 1 -5 cm'^ g~^) 
appropriate to dark clouds from iBianchi et ali (j2003l ). the 
total mass of the cloud is estimated to be ~ 750—1310 Mq. 
This mass estimate can be compared w i th th e gas mass 
of 826 Mq obtained by IWatt fc Mundvl (|1999f ) using the 
C^^O measurements. Assuming 750 M0to be distributed in 
a spherical volume of radius 0.6 pc, we find the mean den- 
sity to be 3.4 X 10"* cm""^. For this mass and mean density. 



the free fall time of the cloud is estimated to be ^ 3 x 10^ 
yrs. 

The visual extinction (in magnitudes) for constant tem- 
perature and dust density distribution along the line-of- 
sight can be obtained by using the following expression: 

where 1^ is the flux density along the line-of-sight, ny 
is the dust emissivity normalised to emissivity in the V- 
band, and Ay is the extinction along the line-of-sight. For 
a dust temperature of 37 ± 5 K at I RAS 18511 core and 
K850fj.m/Kv = 4.0 x 10~^ obtained bv lBianchi et all (|2003( ) 
for dark clouds, the maximum extinction using the peak 
850 /im flux is estimated to be 86 ± 15 magnitudes. The 
average extinction has been estimated assuming that the 
total mass is distributed uniformly in the area enclosed 
within the 10% contour level (^ 1.2 x 1.8 pc^) at 850 /im. 
Assuming a temperature of 20 — 30 K over this region, the 
average value of extinction is estimated to be Ay ^42 — 24 
magnitudes. 

3.2. Mid infrared emission 

The MIPS GAL 70 /im image shows the presence of three 
sources in this region. This image is shown in Fig. O The 
IRAS position is also marked. These sources are better re- 
solved in the MIPS 24 /xm image. For convenience, we call 
these sources A, B and C. The location of these sources 
is shown on the 24 /im image in Fig. [51 While the IRAS 
peak coinicides with A, B lies nearly 20" (0.4 pc) to the 
south-west of A. C is the faintest of the three and is sepa- 
rated by ^90" from A. A is saturated in the MIPS 24 /xm 
image. The fluxes and positions of these sources have been 
extracted using the methods described in the Sect. 2.2.1. 
Table [1] lists the positions as well as fluxes of A, B and C, 
respectively. 

The four GLIMPSE-IRAC images show a bright source 
at A and a small compact group of stars located to the 
south-west of this source (associated with B). Figure [2] 
shows the locations of A, B and C on the grayscale 8.0 /im 
image. The IRAC 3.6 /im image of IRAS 18511 is also shown 
in Fig. [21 The bright source associated with A is saturated 
in all the four IRAC bands. However we have extracted the 
lower limits to the flux of this source using the method of 
aperture photometry as mentioned in Sect. 2.2.2. A compar- 
ison of the MIPS GAL 24 /im image with the higher angular 
resolution IRAC image shows that B consists of a group of 
sources. One can see an extension corresponding to A in 
the 8.0 /im image suggesting that A may consist of two or 
more sources. In the GLIMPSE 3.6 or 4.5 /im images, the 
emission from C is very faint. The emission from C is diffuse 
in appearance and stronger in the 8.0 /im image. 

3.2.1. The filaments around IRAS 18511 

An investigation into the 8.0 /im image shows diffuse emis- 
sion in the neighbourhood of IRAS 18511. It is particu- 
larly interesting to note the filamentary structures (white 
in Fig. [3]) seen in absorption against this diffuse emission. 
IRAS 18511 seems to lie on one clump and extension of 
this clump (in the form of filaments) can be seen towards 
the south-west of A. A clumpy structure can be observed 
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towards the east of IRAS 18511. An estimate of extinction 
towards these filamentary and clumpy structures has been 
obtained by using the following relation: 



I = Ir. 



(3) 



where lo is the background source intensity (determined in 
the vicinity of the filamentary structures), / is the observed 
intensity, and rg is the optical depth at 8.0 yum. One ob- 
serves considerable variations of lo across the image. We 
have therefore assumed a different value of lo for each pixel 
as follows. A median filter image has been constructed us- 
ing a window size of 99 pixels (2'). Using the median value 
at every pixel as the corresponding /q, a visual extinction 
(Ay) map has been constructed taking Ay/rg = 22.5. This 
relation (Ay/rg) has been obtained by taking the ratios 
of extinctions ( Ag/A^ = .43, Ax/Ay = 0.112) obtained 
by Inde betouw et all (|2005l) and lRieke fc Lebofskvl (|l985l) . 
The pixels representing values higher than the median have 
been masked. The Ay values are found to have values in the 
range 6 — 8 in the filaments as well as towards the clumpy 
structures. Converting the Ay map to a column density 
map, the mass in the filaments as well as in the clumps 
have been estimated. While the mass of gas in the clump 
(aj2ooo,i5j2ooo ~ 53"* 47!3 +01°51'0") is found to be 
~ 46 M0, the mass in a few filaments is found to be in 
the range ~ 5 — 15 Mq. The total mass of the filaments 
and clumpy structures shown in Fig. [3] is ~ 120 Mq. These 
mass estimates represent lower limits as Iq has contribu- 
tions from foreground as well as background emission. In 
addition, the total mass estimate also suffers from blank- 
ing of pixels near the 8.0 /im sources. In particular, since 
IRAS 18511 A is very bright at 8.0 fim, the clump contain- 
ing IRAS 18511 A is masked out in the extinction map. 
It is interesting to note that there is an overlap of a fila- 
ment in the 8.0 ^m image with the JCMT-SCUBA 850 /im 
sub-millimetre emission towards the south-west. 



3.2.2. Associated young stellar objects (IRAC Colour-colour 
diagram) 

In order to select IRAC sources from the GLIMPSE cata- 
log for further study, we have defined a 'region of interest' 
around IRAS 18511. We take this to be a region overlap- 
ping the sub-millimetre emission from cold dust. In other 
words, we take the sources in a region enclosed by the 10% 
contour level of the peak of the the sub-millimetre 850 /im 
emission. We find a total of 39 sources (including the sat- 
urated source at A). Of these, six are detected in all four 
IRAC bands. These have been plotted in the IRAC colour- 
colour diagram ([3.6]-[4.5] vs. [5.8]-[8.0]) which is shown in 
Fig. m In the diagram, the solid square approximately de- 
lineates the region occupied by Class II sources whereas 
the dotted square covers the region occupied by the Class 
I models of I Allen et al.l (|2004!) (see their Fig. 4). From the 
colour-colour diagram in Fig. |4l we find five sources lying 
within the boxes representing the region occupied by either 
Class I or Class II sources. We label these sources Gl, G2, 
G3, G4 and G5. Of these 5 objects, G2, G3, G4 and G5 lie 
within the Class I box while Gl lies in the overlap region 
between Class I and Class II. These sources are also shown 
in Fig. [2] and the details of the coordinates and fiuxes of 
these sources are included in Table [21 



3.3. Near infrared emission 

The near infrared Palomar band image of IRAS 18511 
is shown in Fig.[5l Unlike many other star forming regions, 
we do not detect near infrared nebulosity in these images. 
Among the sources extracted from the Palomar images in 
the J, H and bands, we select a sample of sources ly- 
ing within the 'region of interest' (see Sect. 3.2.2 for de- 
tails). We find a total of 68 sources, including the saturated 
sources. The fluxes of three sources which are saturated in 
the Palomar images have been taken from the 2MAS^. 
This includes the source associated with A, which is satu- 
rated in the H and Kg bands of the Palomar images. Of the 
total of 68 sources, we find that 27 sources are detected in 
all three JHK, bands. 



3.3.1. Associated young stellar objects (NIR and IRAC 
colour-colour diagrams) 

In order to identify the young stellar objects in this region, 
we have constructed the colour-colour diagram of J-H vs. 
H-Ks using the 27 sources detected in JHKs bands. This is 
shown in Fig. [5] (left). The loci of the main-sequence and 
giant branches are shown by the solid and dotted lines, re- 
spectively. While t he short-dashed lin e represents the locus 
of T Tauri stars (|Mever et al.l 119971 ). the dot-dashed line 
represents the reddening vector of the main-sequence stars. 
The long-dashed line s hows the locus of Herbig Ae/Be stars 
(jLada fc Adamsl[l992[ ). We have assumed extinction values 
of Aj/Ay=0.282 Aff /A y=0.175 and Aif/Ay=0.112 from 
iRieke fc Lebofekvl (Il98l . The sources having infrared ex- 
cess are those lying to the right of the reddening curve. 
We find a total of 10 sources having infrared-excess. These 
sources are included in Tabled! The source associated with 
IRAS 18511 A has an infrared excess as well as very large 
reddening (J-H~4.1, H-K^ ^2.7) as can be seen in the CC 
diagram. IRAS 18511 A is listed as source number 7 in 
Tabled! 

We have also searched for Spitzer-IRAC counterparts 
to the near infrared objects. We have used a search radius 
of 0"8 . A total of 24 P a lomar sources have IRAC counter- 
parts. IWang fc Loo^ (|2007f) have used J-[3.6] vs. K,-[4.5] 
colour-colour diagram to identify the young stellar objects 
in groups around Herbig Ae/Be stars. In our sample we find 
6 objects that are detected in the J, Kg, IRACl (3.6 /im) 
as well as IRAC2 (4.5 /im) bands. These objects have been 
plotted in the J-[3.6] vs. Ks-[4.5] colour-colour diagram in 
Fig. [5! (right). All the sources in this colour-colour diagram 
have been dereddened by 7 magnitudes of visual extinc- 
tion (1.8 mag kpc~^), corresponding to the extinction due 
to interstellar m e dium (Whittet ,19 92). We have used the 
IWang fc LoonevI (|2007l ) fine (their Eqn. 1) to separate the 
normal stars and the young stellar objects. This is shown 
by the solid line in the Fig. [S! The dotted line represents 
the young stellar object (YSO) locus given by them. From 
this diagram, we find 3 objects which can be characterised 
as YSO candidates based on the above criterion. The de- 
tails of the 3 YSO candidates from the J-[3.6] vs. Ks-[4.5] 
are included in Table d! All the young stellar objects se- 



This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the NASA 
and the NSF. 
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lected from various colour-colour diagrams ([3.6]-[4.5] vs. 
[5.8]-[8.0], J-H vs. H-K^, J-[3.6] vs. K^-[4.5]) are overplot- 
ted on the grayscale Palomar band image in Fig. [5] and 
listed in Tabled 

3.4. Radio emission from ionised gas 

The VLA data show weak emission (point source) at both 
8.5 and 15 GHz in the IRAS 18511 region. The emission 
is from a point source with coordinates (q!,/2000: <5,/20oo 
= ) 18'* 53™ 38^67 +01°50'13'.'0. The location of the ra- 
dio point source coincides with G5, which is marked as 
a cross in the 8.0 /xm image in Fig. O The flux densi- 
ties are 0.69 mJy/beam and 0.68 mJy/beam at 8.5 and 
15 GHz, respectively. The size of the radio emitting re- 
gion is ^5" which corresponds to 0.09 pc at the distance of 
IRAS 18511. The spectral index is -O.OStgJo indicating the 
nature of this emission to be (optically thin) free-free emis- 
sion. Thus, the ionised gas here is from a small extremely 
co mpact region around this source. U sing the f ormulation 
of lSchraml fc Mezgeil (|1969D . as well as lPanagial ()19731 the 
ZAMS spectral type of this source is estimated to be B2- 
Bl (flux of Lyman continuum photons is ~ 1.3 x 10"*^ s~^). 
There is no radio emission detected from the other sources 
in this region, including the bright A source, upto 2a levels 
of 100 /iJy and 300 yuJy at 8.5 and 15 GHz, respectively. 

3.5. Source C 

The source C is clearly detected at 24 and 70 /im and ap- 
pears as a faint diffuse emission. The weak diffuse emis- 
sion shows filamentary structures within it. Also, towards 
the north, we discern a filamentary structure with a sharp 
edge. The cause of this emission is unclear. The emission at 
8.0 /im is stronger than at 24 /im indicating the presence of 
strong Polycyclic Aromatic Hydrocarbons (PAH) emission. 
The angular size of C is ^35" which corresponds to a size of 
0.7 pc at 3.9 kpc (distance to IRAS 18511). While C being 
an irregular galaxy cannot be ruled out, it is possible that 
C is a Galactic object. 

3. 6. Spectral energy distributions 

The fluxes at different wavelengths have been extracted for 
various sources in IRAS 18511 region and these have been 
used in constructing their SEDs. For the sake of compari- 
son, we have also constructed the SED of IRAS 18511 as 
a single source using IRAS as well as MSX fluxes. The lu- 
minosity obtained by integrating the area under the IRAS- 
MSX curve gives an estimate of the total luminosity. Figure 
El (left) shows the SED of IRAS 18511 constructed using 
IRAS and MSX fluxes. Also plotted in this figure is the 
SED of A. It is to be noted that the fluxes of A at the 
IRAC and MIPS 24 /zm band are lower limits as A is sat- 
urated in these wavebands. The luminosity of IRAS 18511, 
obtained by integrating the IRAS-MSX SED, is 1.1 x lO^* 
Lq. The lower limit to the luminosity of A is 7.3 x 10"^ Lq, 
derived using Spitzer photometry and JCMT-SCUBA data. 
Hence, most of the luminosity of IRAS 18511 (at least 66%) 
is due to the protostar(s) in A. 

The SEDs of the individual members (G1-G5) have been 
constructed using wavelengths at which they are resolved 
and detected. These SEDs are shown in Fig. [7] (right) for 



no dereddening applied. The SEDs of these individual ob- 
jects have been constructed using the Palomar as well as 
IRAC fluxes. From Fig. [3 we see that all the SEDs rise 
rapidly with increasing wavelength. In order to compare 
the SEDs relative to each other, we have normalised them 
with respect to the fluxes in H band. The SEDs of Gl, 
G2, G4 and G5 have been dereddened by Ay ~7 magni- 
tudes of visual extin ction (due to ISM ) as well as by Ay 
^ 22 mag derived bv llshii et all (|2002D . Figure [8] shows the 
SEDs relative to each other for a dereddening of Ay ~ 7 
mag (left) and 22 mag (right), respectively. By integrating 
the area under the curves (near to mid infrared), we have 
estimated lower- limits to the luminosities {l-lNIR-8^lm) of 
these objects. These luminosities are listed in Table [3] for 
the sources G1-G5 for dereddening of 7 and 22 mag of vi- 
sual extinction. The source G5 appears to be identical to 
the VLA source discussed in Sect. 3.4 and hence we expect 
it to be of B1-B2 spectral type with an effective tempera- 
ture ^ 20, 000 K and a bolometric luminosity 4400 Lq. 

3.6.1. Robitaille et al. models 

From Fig. [H we see that the sources show a rise in spectral 
energy distributions with an increase in wavelength. Their 
IjNiR-Sfim (Hsted in Table [3]) suggest that at least a few of 
these are massive young objects. In order to get a qualita- 
tive estimate of the evolutionary stage of the cluster mem- 
bers, we have fitted the SED o f few c luster members with 
the models of iRobitaille et all (|2007[ ) (hereafter RWIW). 
They have computed a large set of radiation transfer mod- 
els and obtained the SEDs for a reasonably large param- 
eter space. These SEDs can be fitted to multi-wavelength 
observational data of single sources to constrain the phys- 
ical parameters and the evolutionary stage. However, it is 
important to note that the best constraints are obtained 
if the near, mid and far infrared as well as sub-millimetre 
data are included. 

We have fitted the RWIW models to two sources, G4 
and G5. The models were fit to the SEDs dereddened by 
22 mag of visual extinction. The LAr/ij-g/jmS of G4 and G5 
are the highest among the G sources. The SEDs have been 
constructed using near infrared Palomar and mid infrared 
IRAC data. Upper limits at MIPS 24 /im (corresponding 
to the IRAS 18511 B) and at JCMT-SCUBA 850 /xm (flux 
density at the position of IRAS 18511 B) have also been 
used as constraints to the modelling. It should be noted 
that our sub-millimetre angular resolution is not sufficient 
for this purpose and hence the fits are non-unique. We have 
also put additional constraints based on luminosity and 
mass of envelope/disk (85 Mq based on the sub-millimetre 
850 /im flux at the position of IRAS 18511 B). The models 
are selected based on the least-squares chi-square method. 
For G5, additional constraints are available since it is esti- 
mated to be of ZAMS spectral type B1-B2 (based on the 
VLA observations). Among the four models for G5, Model 1 
simulates a young protostar and the Te/ / of the central ob- 
ject does not produce Lyman continuum photons to ionise 
the surrounding gas. On the other hand, the other models 
assume more evolved central stars (Te// ^20000 K). Figure 
[9] shows some sample RWIW models for G4 and G5 along 
with the observed SEDs. The parameters are listed in Table 
m In the Table, column 1 lists the source (G4/G5) column 
2 (M*) lists the mass of central object, column 3 lists the 
effective temperature of the central objects, column 4 lists 
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the total luminosity, column 5 lists the inclination angle 
with respect to the observer (90°is edge-on), column 6 (M) 
represents the envelope accretion rate, column 7 (Me„^) 
lists the mass of envelope, column 8 (M^isfe) hsts the mass 
of disk , column 9 (Ay) lists the extinction along the line- 
of-sight and column 10 lists the age. The luminosity listed 
in the Table is the total (non-isotropic) luminosity (central 
object + envelope and/or disk) computed by RWIW. 

It is interesting, especially for G5, that the observed 
fluxes in the mid infrared (IRAC) bands have a steeper 
wavelength dependence than predicted by these models 
which perhaps underestimate the extinction. Another likely 
possibility is that emission from small grains or from 
Polycyclic Aromatic Hydrocarbons (PAH) effects the obser- 
vations. The results indicate that the SEDs of few sources 
can be explained by massive Class I and Class H type 
sources i.e., sources with the massive central objects (7— 10 
Mq) alongwith envelopes/disks. Thus, the models of low- 
mass objects extended to more massive candidates can rea- 
sonably explain the SEDs of some cluster members of IRAS 
18511. 

4. Cluster simulations 

In the IRAS 18511 region, a number of sources have been 
detected in the near and mid infrared. In regions of high 
extinction like IRAS 18511, these sources could be highly 
reddened luminous objects or low-luminosity low-extinction 
objects. The objective of the cluster simulations is to inves- 
tigate this. We simulate a cluster of young objects embed- 
ded in a cloud of gas and ascertain whether the observations 
of IRAS 18511 in various near and mid infrared bands are 
in qualitative agreement with the predictions of the model 
as our statistics are too low for a quantitative comparison. 
We explore the cluster membership in terms of the frac- 
tion of objects of different evolutionary stages (Class I and 
Class II). As we are exploring a young embedded cluster, 
we have not considered objects in the Class III phase or 
later. This modelling will allow us to obtain a qualitative 
estimate of the evolutionary stage of the cluster by varying 
the fractions of Class I and Class II sources in the cluster. 

The inputs to these simulations include (a) observables 
- quantities that have been incorporated in the model based 
on observations, and (b) assumptions. The observables in- 
clude (1) limit on the most massive object in the cluster, 
(2) the total luminosity of the cluster, (3) mass and size of 
the spherical cloud of gas in which the cluster is embed- 
ded, and (4) the distance. We have assumed the following: 
(1) Initial Mass Function (IMF) and (2) Star formation his- 
tory (SFH) . The output that can be compared directly with 
the observations is the number of detected cluster members 
in the Ks (2.12 ^m), IRACl (3.6 fim), IRAC2 (4.5 fim), 
IRAC3 (5.8 Atm), IRAC4 (8.0 ^lul) and MIPSl (24 /im) 
bands for given completeness limits. A detailed description 
of the model is given in Appendix B. 

4.1. Results 
4.1.1. Model 

The upper mass limit of the cluster has been taken to be 
10 Mq corresponding to the ZAMS spectral type of a sin- 
gle object inferred from the bolometric luminosity of IRAS 
18511. The model has been run for the Salpeter IMF. The 



mass and radius of the spherical cloud of gas have been 
taken to be 750 Mq and 0.6 pc, respectively. These have 
been derived from the 850 /im JCMT-SCUBA observations. 
This corresponds to a maximum visual extinction of 69 
magnitudes (including the Ay ~ 7 magnitudes due to the 
ISM). The total luminosity has been taken to be 1.1 x 10'* 
Lq which is the total luminosity of IRAS 18511. The model 
was run 1000 times for each of the following cases of star 
formation history: 

1. Coeval formation with an age of 0.5 Myr correspond- 
ing to Class I type sources - Every cluster member is 
assumed to have a Class I type spectrum. 

2. Coeval formation with an age of 1 Myr corresponding to 
Class II type sources - Every cluster member is assumed 
to have a Class II type spectrum. 

3. Uniform star formation rate between 0.01 — 1 Myr - This 
corresponds to the general case incorporating Class I as 
well as Class II objects. For every cluster member of a 
certain mass, an age is randomly assigned. If this is less 
than 0.5 Myr, its spectrum is taken to be of the Class I 
type. For ages greater than 0.5 Myr, the cluster member 
is a Class II type object. 

The median number of cluster members is ~ 300 while 
the median mass of the cluster is 230 Mq. We have also 
estimated the magnitude distribution from the cluster for 
the case of a uniform star formation rate (Case 3). This is 
shown for Kg band apparent fluxes in Fig.[Tn]for the general 
case incorporating Class I and Class II objects with ages 
between 0.01 and 1 Myr. 

4.1.2. Comparison with observations 

A comparison of the model results with the observations 
has been carried out in terms of magnitude distribution 
(Fig. (TU]) and number of sources detected (Fig. [Tl|). Figure 
[TO] shows the simulated (Case 3) as well as observed mag- 
nitude distributions from the Palomar NIR data. The ob- 
served sources used in the figure are the young stellar ob- 
jects detected in the band and listed in Table [H Also 
marked on the figure is the Palomar band sensitivity 
limit. A comparison of the observed and simulated magni- 
tudes shows that the sources detected in observations are 
brighter than those simulated. This could be because the 
density distribution of the cloud is inhomogeneous while 
the cluster simulations assume a constant density distribu- 
tion of gas. Another possibility is that the distribution of 
young stellar objects is more extended than the size of the 
cloud. From the simulations, it is clear that with the sensi- 
tivity limit of Palomar, only a small fraction of sources in 
the cluster are detected. 

From the model, the median number of detected sources 
in each band has been derived based on the completeness 
limit of the instrument used for comparing the results. 
Figure [TT] shows the number of detected cluster members 
predicted by the model in each band. The errorbars indicate 
the quartile values encompassing 50% of the number of de- 
tected cluster members. The cross and dotted line represent 
Case 1 (Class I); open circle and dashed line represent Case 
2 (Class II); and, filled square and dot-dashed hue represent 
Case 3 (general case incorporating Class I Class II objects). 
The observed number of sources (from Palomar image and 
IRACl, IRAC2, IRAC3 & IRAC4 bands of GLIMPSE im- 
ages as well as from MIPS 24 /im of MIPSGAL image) are 
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shown by the sohd circles and sohd lines. The wavelengths 
for different cases have been slightly shifted for better view- 
ing. 

The observed number of sources in each band is shown 
by a solid line depicting a range of possible number of clus- 
ter members, between the lower and upper limits. The lower 
limit is obtained from the number of young stellar objects 
detected in that band within the 'region of interest' by using 
the colour-colour diagrams (ref Table [2]). The upper edge 
of each solid line represents the total number of sources de- 
tected within that band. For the case of the MIPS 24 /im 
band, there are only two sources detected within the region 
of interest. IRAS 18511 B shows an extension at 24 /zm and 
therefore we have plotted the number of sources detected 
as 2. Obviously, this represents a lower limit as the sources 
are unresolved due to lower angular resolution. 

From Fig. [Til we find that among all the three cases, 
the general case incorporating Class I ('^50%) and Class 
II (~50%) sources agrees well with the observations. The 
models with Class I and Class II sources alone are inconsis- 
tent with the observational data. We, therefore infer that 
the sources in IRAS 18511 are very young and deeply em- 
bedded. In Ks as well as in IRACl bands, the general case 
(Case 3) as well as Class II (Case 2) model results are con- 
sistent with observations. This is expected as more Class 
II sources should be detected in the band. However, in 
the IRAC2 and IRAC3 bands, all the model results agree 
with observations within the quartile values. For the IRAC4 
band, the Class I as well as the general case agree with ob- 
servations better than with the Class II case. Finally, at 
24 fiui, the number of detected sources from the Class I 
model and the general case are higher than that detected 
from the Class II model. It must be noted that the obser- 
vations at 24 /im represent a lower limit due to the low 
angular resolution of the MIPS instrument (^6"). Among 
all the cases, the general case (Case 3) assuming a mix of 
^50% Class I and ^50% Class II sources can be said to be 
a reasonable fit to observations. The model suffers from a 
few limitations which are listed below. 

— The cluster is embedded in a cloud of gas which is homo- 
geneous and clumpiness/density distributions has not 
been taken into consideration. 

— The number of detected cluster members from the sim- 
ulations (which are point-like) is compared with the ob- 
served number of sources. A better comparison would be 
possible if the results of the simulations were convolved 
with the beam size of the instrument. 



5. Discussion 

5.1. Clustering 

We have attempted to study the clustering of the sources 
around IRAS 18511 using the NIR Palomar im a ges. T he 
method used is described in detail bv lTesti et al.l ( 19981 ) to 
study the clustering around Herbig Ae/Be stars. A density 
profile is computed centred on the bright IRAS 18511 A 
source to look for evidence of clustering. Density profiles 
are computed using the Palomar Ks, IRAC2 (4.5 ^m) and 
IRAC4 (8 /im) band sources as shown in Fig. [T^l In this 
plot, the star densities derived for the IRAC2 and IRAC4 
bands have been scaled by factors of 1.7 and 6.4, respec- 
tively for comparison with the radial profile derived for the 



Ks band. Further, the radii of the annuli of the IRAC bands 
are slightly shifted for better viewing of the plot. From the 
profile, we do not see any 'clear' evidence for clustering cen- 
tred on A. This is consistent with the results of the previous 
sections, i.e. the extinction due to the associated molecu- 
lar cloud is rather high and there may be other embedded 
sources. This is consistent with the results of our simula- 
tions, i.e. only a small fraction of sources are detected. 

5.2. Extinction at 8.0 fim vs. 850 fim 

A comparison of the morphology of the sub- millimetre dust 
emission and the filamentary structures seen in the 8.0 /im 
emission map indicates that these filamentary structures 
comprise of cold dust. Figure [3] shows the 8.0 /im grayscale 
image with 850 /im contours overlaid. The morphology of 
the 850 /im emission traces the filamentary structures close 
to IRAS 18511. It is evident that these filaments represent 
the cold dust comprising the infrared dark cloud with IRAS 
18511 located at the peak of one clump. 

We have compared the extinction values obtained from 
the JCMT-SCUBA 850 /tm image (using Eqn. 2) with that 
from IRAC 8.0 /im image (using Eqn. 3). Resampling the 
images and a pixel-to-pixel comparison of the extinction 
values in the region of overlap of the filament (to the south- 
west of IRAS 18511 A) shows that the visual extinction 
values obtained using the JCMT-SCUBA sub-millimetre 
image is ~ 8 times the the Ay values obtained from the 
IRAC 8.0 /(m image. However the Ay values quoted for the 
IRAC 8.0 /tm image are lower limits for the following rea- 
son. Equation 3 assumes the emitting dust to be behind the 
absorbing filaments and this is clearly not correct since lo 
has contribution from foreground emission. Consequently, 
the values of Ay derived are an underestimate. And also, 
the mass estimates derived fo r the filament s and the clump 
in Sect. 3.2.1 are lower-limits. i Brand et al.l (j200lD have im- 
aged a number of molecular lines in the IRAS 18511 region 
(i^CO, CS, HCO+). Of particular interest is the HCO+ 
emission integrated over the red part of the emission. This 
emission extends along the direction of the filaments to the 
south-west, as seen in the JCMT-SCUBA 850 /im image. 
Further, the emission is at a velocity (V/sr) of ~ 59 km s~^. 
This is the V/sr of IRAS 18511 and hence, the filamentary 
structures are associated with the molecular cloud of IRAS 
18511. 

Using the fact that the Ay values derived from the 
850 /tm sub-millimetre map is '^8 times larger than that 
obtained from the IRAC 8.0 /tm map, we can estimate the 
fraction of foreground emission contributing to the total 
emission detected (in the neighbourhood of the filaments). 
We estimate that this fraction of foreground emission is 
^^80% of the detected emission. In other words, most of 
the emission is foreground emission. This is not surprising 
considering that IRAS 18511 is a distant object (~3.9 kpc). 

5.3. Evolutionary stage of cluster members 

It is interesting to compare the multiwavelength observa- 
tions of A and B, although each comprises of more than one 
source. A is the brightest source in the near, mid and far in- 
frared as well as in the sub-millimetre wavebands. However, 
there is no detected radio emission associated with A. On 
the other hand, G5 (associated with B) has free-free ra- 
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dio emission associated with it. The peak of molecular Hne 
(C^^O and ^^CO) emission, as observed bv lWatt fc Mundvl 
(|l999t ) is located near B and one can see an extension cor- 
responding to A. The absence of radio emission near A and 
other observations suggest the following possibilities: (1) 
A is in an earlier evolutionary stage than B and the ultra- 
compact Hii region has not yet formed. This would explain 
the absence of radio emission as well as their brightness at 
infrared wavebands (2) A being more massive than B, its 
ultracompact Hii region i s being choked by the infalling 
matter through accretion (|Ketoll2003D . iKurtz et all (|2004 ) 
have carried out a VLA survey in the 44 GHz (Class I) 
methanol maser line in several star forming regions includ- 
ing IRAS 18511. They find a maser in the vicinity of IRAS 
18511 A (synthesized beam size ^ 2'.'2 x 1'.'5). The posi- 
tion of the maser (positional uncertainy ~ 0'.'5) is shown in 
Fig. [21 They also find in their survey that Class I methanol 
masers are typically offset by ~ 0.2 pc (median value) from 
other massive star formation signposts like Hii regions or 
water masers. ICveanowski et all (|2007l ) have studied the 
massive protocluster associated with S255N and they find 
that the Class I methanol are aligned in the direction of the 
outfiow. In the case of IRAS 18511, the methanol maser is 
at distance of 0.1 pc from IRAS 18511 A and 0.4 pc from 
G5 (VLA source). Also, no water masers or H ii region have 
been detected near IRAS 18511 A. Further, it is believed 
that shocked gas from outflows /accre tion might be respon- 
sible for Class I methano l masers JP lambec k fc MentenI 
119901: iJohnston et al.|[T99l ). IZhang et al. (200^3) have im- 
aged IRAS 18511 (-30") in the CO (J=2 ^ 1) line and 
the emission shows a peak close to IRAS 18511 A. This 
is unlike the case of C^^ O where the peak of e mission is 
located at IRAS 18511 B (|Watt fc Mundvl 1 19991) . There is 
also evidence for an outflow from IRAS 18511 A since wing s 
are seen in the line proflles observed bv lBrand et al.l (|200lh . 
All these facts suggest that IRAS 18511 A is in an earlier 
evolutionary state than IRAS 18511 B. 

Considering the above observations, analyses and simu- 
lations of the cluster associated with IRAS 18511, the fol- 
lowing likely scenario emerges. This is a young cluster with 
sources in an early evolutionary phase (Class I and Class 
II). The large extinction indicates that probably there are 
more low-mass cluster members hidden in the molecular 
cloud which will start appearing as the cluster evolves and 
the molecular cloud disperses. IRAS 18511 A is likely to 
be an intermediate mass/massive young object, probably a 
Class I object. Its luminosity indicates that it is probably 
the most massive object in the cluster. On the other hand, 
the presence of radio emission near G5 indicates that G5 
could be an older member of the cluster. The absence of 
other young objects in the vicinity of IRAS 18511 A can 
possibly be explained by large extinction as indicated by 
the JCMT sub- millimetre maps (peak of emission is at A) . 
Further study of IRAS 18511 (high angular resolution mid 
infrared maps and molecular line maps) can shed light on 
the likely scenario. 



5.4. Precursor to Herbig star type cluster? 

The clusters associated wit h Herbig stars and their prop- 
erties h ave been studied by T esti et all ()1997f ). iTesti et ahl 
(|1998f) . iTtesti et al . (1999). In particular, the mid infrared 
emission (10 /xm) from a few such clusters has been investi- 
gated bv lHabart et al.l ()2003f ). The luminosity as well as the 



clustering suggests that IRAS 1 8511 could be a precursor 
to a Herbig cluster. iTesti et all (|1999 f) report the tendency 
of early Herbig Be stars to be surrounded by dense clusters 
of lower mass companions. Some of the clusters associated 
with earlier Herbig stars (of ZAMS spectral type BO) are 
found to have 75 members, where the lower mass limit is 
~ 0.2 Mq for 2 magnitudes of extinction in K. It is interest- 
ing to compare these numbers with the model results which 
suggest cluster membership to be ~ 200 for masses greater 
than 0.2 Mq. While the model results predict the number 
of cluster members to be — 2 — 3 times the number observed 
b v ITesti et a l.' (1999), it is useful to note that the clusters 



of ITesti et al.. (1999) are more evolved as compared to the 
model which is an embedded cluster. And evolved clusters 
with intermediate m ass Herbig Ae/Be sta rs retain less than 
50% of its members (IWeidner et al.|[2007l ) Further, the size 
of the embedded cluster (size of cloud from sub-millimetre 
map) associated with IRAS 18511 is —0.6 pc. This com- 
pares well with the cluster sizes of 0.2-0.7 pc obtained for 
clusters associated with Herbig Ae/Be stars of spectral type 
BO (Testi et al. 1999). The mid infrared (excess) emission 
as well as the large extinction (it is not visible optically) 
indicates that the main source (A ) is in an early evolu- 
tionary stage. iHabart et al.l (|2003l ) flnd objects with mid 
infrared emission from circumstellar disks and envelopes in 
flve out of twelve flelds of known Herbig AeBe stars they 
studied. The IRAS 18511 field also shows many objects in 
mid infrared. In short, the following considerations: (a) lu- 
minosity (b) evolutionary stage of the objects (c) size of 
associated cluster and (d) number of cluster members, sug- 
gest that the source associated with IRAS 18511 is a proto- 
cluster associated with a candidate precursor to a Herbig- 
star cluster . This is i n accordance with the suggestions of 
ITesti et al.l (|1998[) and lMolinari et al.l (|2000l ). We therefore, 
believe that IRAS 18511 represents an early stage of clus- 
ters associated with an intermediate mass object. 



6. Summary 

With the aim of studying the early evolution of clusters, we 
have selected IRAS 18511 to carry out a detailed investi- 
gation. Using emission at sub-millimetre (JCMT-SCUBA), 
infrared (Spitzer-MIPS, Spitzer-IRAC and Palomar) and 
radio (VLA) wavelengths, we have studied the main source 
as well as the cluster associated with IRAS 18511. The re- 
sults are consistent with simulations of a young embedded 
cluster incorporating Class I and Class II sources. Based 
on the luminosity, and properties of the associated cluster 
(number of members, size of cluster and evolutionary stages 
of objects), we conclude that IRAS 18511 is a protocluster 
with the most massive object being a candidate precursor 
to a Herbig Ae/Be star. 
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Table 1. The MIPS flux density details of sources in IRAS 18511. 



Source 


Position 


Flux Density (Jy) 




J2000 


24 


70 /im 


A 


18'' 


53"' 


37!85 +01" 


50' 


30'.' 5 


> 20.6 


145.7 


B 


18'' 


53"' 


38^44 +01^* 


50' 


15'.' 1 


6.4 


61.2 


C 


18'' 


53"' 


4S55 +01^* 


49' 


50'.'4 


0.8 


22 



Table 2. Young stellar objects around IRAS 18511 selected from Palomar (NIR) as well as Spitzer-IRAC (MIR) colour-colour diagrams. 
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G034.8194±00.3420 


13.42 ±0.17 
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283.415892 


1 


.834844 


18, 


,02 
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0, 


,16 


16.28 ± 0.06 


14.66 ± 0.04 


G034.8187+00.3417 


13.15 ±0.14 


12.81 ±0.12 







" Source saturated in tlie GLIMPSE-IRAC images 
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Table 3. Luminosities of the sources G1-G5 in IRAS 18511 region detected in all four IRAC bands 



S. No. 


Luminosity {hNiR-s^im) 




Lq {Av ~ 7 mag) 


Lq (Av ~ 22 mag) 


Gl 


7 


16 


G2 


2 


9 


G3 


19 


36 


G4 


184 


412 


G5 


238 


492 



Table 4. Parameters of the iRobitaille et all ()2007D models for G4 and G5, plotted in Fig. [H Column 2 (M,) Hsts the 
mass of central object, column 3 lists the effective temperature of the central objects, column 4 lists the total luminosity, 
column 5 lists the inclination angle with respect to the observer (90°is edge-on), column 6 (M) represents the envelope 
accretion rate, column 7 (Me„„) lists the mass of envelope, column 8 (Mdisk) lists the mass of disk , column 9 (Ay) lists 



the envelope extinction along 


; the line-of-sij 


;ht and column 10 lists the age. 










Obj. 


M. 


Te// 


Luminosity 


Incl. angle 


M 




Md.sk 




Av 


Age 




(Mo) 


(K) 


(Le) 


(deg) 


(Me/yr) 


(Mo) 


(Mo) 




(mag) 


(yr) 


G4 Model 1 


6.7 


4241 


5.2 X 10^ 


18 


3.9 X lO"'' 


2.4 


1.6 X 10" 


-'2 


19.3 


5 X 10'' 


G4 Model 2 


9.2 


24552 


4.6 X 10^ 


81 







2.1 X 10" 


-2 


24 


1 X 10® 


G5 Model 1 


10.1 


4351 


2.2 X 10^ 


18 


1.1 X 10"^ 


47.5 


4.0 X 10" 


-2 


1960 


4 X 10^ 


G5 Model 2 


7.5 


21770 


2.3 X 10^ 


87 


1.1 X 10"^ 


1.3 


8.4 X 10" 


-3 


414 


4 X 10^ 


G5 Model 3 


7.3 


21262 


2.0 X 10^ 


87 


8.6 X 10"^ 


0.07 


2.9 X 10" 


-5 


12 


1 X 10® 


G5 Model 4 


8.1 


22602 


2.9 X 10^ 


87 







1.6 X 10" 


-3 


5.5 X 10"* 


4 X 10® 



The cavity angle in Model 1 for both G4 and G5 is 2° - 6° while the cavity angle in Models 2 and 3 of G5 is 40° - 50°. 
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RA (J2000) RA (J2000) 

Fig. 1. (Left) JCMT-SCUBA 450 emission from IRAS 18511 region. The contours levels are from 2.0 to 13.4 Jy/beam 
(peak flnx density ^ 13.4 Jy/bcam) in steps of 1 Jy/beam. where the beam size is 10". (Right) JCMT-SCUBA 850 /im 
emission from IRAS 18511 region. The contours levels are from 0.18 to 2.28 Jy/beam (peak flux density ~ 2.28 Jy/beam) 
in steps of 0.2 Jy/beam, beam is 15'.'5. 
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Fig. 2. (Top left)MIPS 70 /xm image of IRAS 18511 region. The plus symbol and the ellipse represent the IRAS position 
and the corresponding l-cr position error ellipse. (Top right) MIPS 24 /xm image of IRAS 18511 region with the sources 

A, B and C marked. (Bottom left) IRAC 8.0 /im image of IRAS 18511 region. The white cross marks the position of 
radio emission at 8.5 and 15 GHz. The solid triangle marks the position of the methanol maser. (Bottom right) IRAC 
3.6 yum image with sources G1-G5 marked. 
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1S''53"'50' 



RA (J2000) 



Fig. 3. Grayscale image of IRAS 18511 at GLIMPSE- IRAC 8.0 /im showing the dark filaments against diffuse emission. 
Overlaid as contours is the JCMT-SCUBA sub-millimetre emission at 850 fim. The region covered by the SCUBA map 
is also shown. 



• Observed 

o Simulations 

(witiiout extinction) 



in 



o Class 1 



CD 

CO 1 
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•Gl 

o Clas: 



C3 



[5.8] - [8.0] 



Fig. 4. Spitzer colour-colour diagram of the observed sources and the models used in the simulations. The filled circles 
represent the sources detected in all four IRAC bands. The empty (red) circles represent the colours of the unreddened 
Class I and Class II models used in the simulations (see text for details). The solid line square approximately delineates 
the region occupied by Class II sources whereas the dotted-line square covers the region occupied by the Class I models 
shown in the models of Allen et al. (2004). The arrow represents the extinction vector for Ay =20 magnitudes. 
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Fig. 5. Ks band Palomar image of the region around IRAS 18511. The axes are in J2000 coordinates. The sources marked 
by open (red) circles are the YSO candidates selected from the various colour-colour diagrams (see Tabled]). The (blue) 
solid triangles are the sources G1-G5. 




Fig. 6. Colour-colour diagrams (CCDs) of sources detected in Palomar and Spitzer-IRAC images within 10% contour 
level of the sub-millimetre 850 /im peak at IRAS 18511 A (see text for details). (Left) J-H vs. H-K^ CCD where the 
locii of the main sequence and giant branches are shown by the solid and dotted lines respectively. The short-dash line 
represents the locus of T-Tauri stars. The dash-dotted straight line follows the reddening vectors of main sequence stars 
(or dwarfs). The long dashed line represents the locus of Herbig Ae/Be stars. The position of IRAS 18511 A is marked as 
'A' in the CCD. (Right) J-[3.6] vs. Ks-[4.5] CCD where the solid line repr esents the line separatiii g the normal stars and 
the young stellar objects while the dotted line represents the YSO locus of IWang &: LoonevI (|2007t ). The arrow represents 
the extinction vector Ay =5 magnitudes. 
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Fig. 7. (Left) The solid line represents the spectral energy distribution (SED) of IRAS 18511 constructed using the 
IRAS and MSX fluxes. The dotted line depicts the lower limit to the spectral energy distribution of A constructed using 
Palomar, Spitzer-IRAC, Spitzer-MIPS and JCMT-SCUBA. While the filled circles represent the fluxes of IRAS 18511 
A, the solid triangles represent the lower limits. (Right) The SEDs of young stellar objects, G1-G5, detected using IRAC 
colour-colour diagram. 




Wavelength (^im) Wavelength (/^.m) 



Fig. 8. Normalised spectral energy distribution (with respect to H band) of sources Gl, G2, G4 and G5, dereddened 
by Av ~7 magnitudes (left) and Ay ~22 magnitudes (right). The solid lines represent the normalised spectral energy 
distributions of the Class I (red) and Class II (blue) sources adopted in the cluster simulations. 
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0.01 0.1 1 10 100 1000 10" 0.01 0.1 1 10 100 1000 10" 
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Fig. 9. Spectral energy distributions of G4 (left) and G5 (right) with some sample models of iRobitaille et all (|2007| ) 
(RWIW). The observed fluxes are shown by filled circles, the triangles denote upper limits to the observed fluxes and 
the lines represent the RWIW models. The parameters of these models are listed in Table ID 
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Fig. 10. Plot of magnitude distribution in the Kg band from the simulations and observations (Palomar). The solid line 
represents the magnitudes of YSOs within the 'region of interest' (Table [2]). The dotted line represents the apparent 
magnitudes from one run of model (general case of constant star formation rate with ages between 0.01 and 1 Myr) 
simulations. The dashed line denotes the Palomar sensitivity limit. 
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Fig. 11. Plot of the detected number of sources as a function of wavelength for cluster model simulations with observa- 
tions. The filled circles and solid lines represent the observations and error bars, respectively. The observation at 24 /im 
represents a lower-limit. The solid squares represent median values of the general case of simulations incorporating Class 
I as well as Class II objects for the general case of constant star formation rate with ages between 0.01 and 1 Myr and the 
dot-dashed line represent the quartile values on either side of the median. The cross symbols and dotted line represent 
the Class I model simulations for an age of 0.5 Myr. The open circles and dashed line denote Class II model simulations 
for 1 Myr. The wavelengths for different cases have been slightly shifted for better viewing. Further details can be found 
in the text. 
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Fig. 12. The clustering indicator (radial profile) centred on the bright NIR source in IRAS 18511, constructed using the 
Palomar (solid circles), IRAC 4.5 /im (solid triangles) and IRAC 8.0 /im (solid squares) band fluxes of the sources in 
the region around IRAS 18511. The star densities derived for the IRAC2 and IRAC4 bands have been scaled by factors 
of 1.7 and 6.4, respectively for suitable comparison with the radial profile in the band. The radii for different cases 
have been slightly shifted for better viewing. 
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Appendix A: Photometry of Spitzer images 

IRAS 18511 consists of a very bright source and a small 
group of stars near it. We, therefore, first explored our pho- 
tometry techniques on a field consisting of isolated sources 
(which are not saturated) with good signal-to-noise. The 
following methods were used. 

1. APEX - The Spitzer Astronomical Point Source 
Extraction (APEX) is a software provided by the 
Spitzer Science Centre for point source extraction and 
photometry. We have used the Linux version MOPEX 
030106 for our study. In this software, the sources are 
first detected using non-linear matched filtering and im- 
age segmentation. Subsequent point source extraction 
is performed by fitting the Poi nt Response Functio n 
(PRE) to the detected sources (jMakovoz et al.ll2004 . 
More details can be found in the APEX Manual. For 
PRE fitting, the PRE has been taken from the Spitzer 
Science Center (SSC) webpage. In addition, aperture 
photometry can also be carried out. We have compared 
the fluxes from apertures of radii 6", 9" and 12". 

2. DAOPHOT - We have used the IDL version of the 
IRAE-DAOPHOT package in order to carry out source 
extraction and aperture photometry with an aper- 
ture radius of 6" and background annulus of 12". 
Appropriate aperture corrections have been applied. 

3. SExtractor - Source Extractor is a program which can 
be used to extract sources from moderately crowded 
fields. The background in this package is, however, de- 
termined from the full image and the objects above a 
given threshold are detected and de-blended. Eor pho- 
tometry, we have used an aperture radius of 6". 

4. AIPS/GILDAS - We have also used the AIPS/GILDAS 
software to compute the integrated fiux down to 10% 
of the peak value after subtracting the background. 
Appropriate aperture corrections have been applied de- 
pending on the size of the integration area. 

All the above methods did a reasonably good task of 
source extraction and photometry. The values of fluxes of 
the isolated bright sources obtained from the methods men- 
tioned above were compared. The values were found to be 
within 5% of each other. We then proceeded to extract and 
carry out photometry of the sources in the region around 
IRAS 18511 which has bright saturated sources located 
close to each other as well as some diffuse emission. The 
saturated pixels were replaced by an average value of the 
neighbouring non-saturated pixels. For this region, we find 
the following: 

- The source extraction by SEXTRACTOR is better 
than that by the PRE fitting methods, APEX and 
DAOPHOT. 

- The aperture photometry methods yield consistent re- 
sults within 5% of the flux values. 

- The PRE fitting method of APEX finds quite a large 
number of sources near the bright source and the flux 
is distributed among them. 

- The flux values obtained by AIPS/GILDAS which inte- 
grates fluxes in circular areas down to 10% of the peak 
brightness after subtracting the background gives con- 
sistent results with those from aperture photometry. 



Appendix B: Description of the cluster simulations 

For a given IMF (Scalo or Salpeter or Kroupa), the clus- 
ter members are randomly selected by using a Monte-Carlo 
method. The lower mass limit is taken to be 0.1 Mq. An 
age (star formation hisory) is assigned to each cluster mem- 
ber according to either of the two prescriptions: (a) coeval 
formation - every cluster member has the same age, or (b) 
uniform formation rate - an age is assigned to every cluster 
member which falls randomly between the ages representing 
the start and finish of the star formation process. The lu- 
minosity and effective temperature of every cluster member 
having a certain mass an d age is determin e d usin g the pre- 
mam sequence tracks of iPalla fc StahleH (|1999D . Because 
of uncertainties in the early pre-main sequence evolution 
of stars, we have assumed a lower limit on the age. The 
pre main sequence track corresponding to 0.5 Myr is ap- 
plied even to objects with ages less than 0.5 Myr. For ages 
of a given mass (as well as masses for a given age) greater 
than that for the pre-main sequence tracks, the luminosities 
and effective temperatures have been obtained by assuming 
them to be lying on the tracks of zero-age main sequence 
(ZAMS) stars. For such Class II objects, it is assumed that 
although the central source has reached the ZAMS stage, 
there is a remnant disk around it. The cluster members 
are randomly added to the cluster till the total luminosity 
reaches the bolometric luminosity of the cluster. It is impor- 
tant to note that the the luminosity of each cluster member 
includes luminosities from the central object (photosphere 
for the Class II objects) as well as from the circumstellar 
material (disk for Class II and envelope for Class I objects). 

The cluster members are now placed at random loca- 
tions in a homogeneous spherical cloud of gas of assumed 
mass and size (derived from the sub-millimetre/millimetre 
maps). In other words, the distribution of the cluster mem- 
bers follows the gas distribution. Thus, there is a column 
of gas in front of every embedded cluster member which 
determines its extinction. Therefore, different cluster mem- 
bers suffer different amounts of extinctions based on their 
positions within the cloud of gas. An additional visual ex- 
tinction of 7 magnitudes (due to ISM) has been added to the 
extinctions. Taking the total extinction into consideration, 
the apparent magnitude of every cluster member is deter- 
mined in each band and is considered detected if the object 
is brighter than the sensitivity limit in that band. The sen- 
sitivities in various bands as well as the extinction laws 
used are listed in Table IB. 11 The near infrare d extin ction 
laws have been taken from iRieke fc Lebofskvl (|l985f ). For 
the Spitzer-IRAC ban ds, we have used the ave rage extinc- 
tion value obtained bv llndebetouw et"an (|2005[ ). At 24 ^m, 
we have used an extrapolated extinction value. 

Every cluster member is classified as Class I or Class 
II. Sources younger than 0.5 Myr are considered as Class 
I and those older than 0.5 Myr are considered to be of 
Class II type. We have taken the age of 0.5 Myr between 
the Class I and more evolved objects as the pre-main se- 
quence ev olution models show uncertainities below the age 
of 1 Myr (|Baraffe et al.ll2002D . For Class I sources, the fiux 
in various bands has been determined based on the spec- 
trum of the well-studied Class I type low mass protostar, 
L1551-IRS5. ICohen fc Schwartd ()1983f ) list the magnitudes 
of L1551-IRS5 in 11 bands ranging from 1.2 /im to 19 fim. 
The ratio of the absolute fiux (corrected for distance of 
L1551) to the luminosity of L1551-IRS5 (^ 35 L©) has 
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been used as a scaling factor to obtain the fluxes of Class 

I type objects in the simulations corresponding to the lu- 
minosities of these objects. The ratios of fluxes (in various 
bands) to the total luminosity are tabulated in Table IrTI 

In order to determine the emission from pre-main se- 
quence objects in the Class II phase, we have used the me- 
dian spectrum of Class II type objects in th e Taurus star 
forming region (Table 6 of lFurlan et"aIll2006D . This median 
spectrum has been constructed using fluxes from bands in- 
cluding the near infrared (JHK^) and Spitzer-IRAC (3.6, 
4.5, 5.8 and 8.0 /im) bands. This median spectrum implies 
flux ratios (flux in various bands with respect to the flux in 
the H-band) which are given in Table EH The normaliza- 
tion has been carried out with respect to the flux values in 
H (1.65 fim) band sinc e the H-band flux is photospheric for 
most Class II objects (jFurlan et al.ll2006l ). As the Class I 
and Class II spectra are of low-mass objects, we have used 
these in the model simulations for objects of all masses. 
Although the pre-main sequence evolution of massive stars 
is not clearly understood, we use these Class I and Class 

II spectra since fits from RWIW models (Sect. 3.6.1) in- 
dicate that the SEDs of cluster members of IRAS 18511 
are well-flt by models comprising massive central objects. 
The Class I and Class II type spectral energy distributions, 
used in the model, are shown in Fig [51 These spectra have 
been normalised with respect to the H (1.65 /im) band. The 
Spitzer colours of Class I and Class II SEDs considered for 
modelling are shown in Figd) 
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Table B.l. Extinction laws and various other constants used in the simulations. F\/Fh represents the ratio of flux 
in each band with respect to t he flux in H b and assumed for Class II sources in the simulations based on the median 
spectrum of Cl ass II sources by Furlan et al.l (p006) and A/L represents the ratio of flux in each band to the luminosity 
of L1551-IRS5 (|Cohen fc Schwartzlll983D " used to estimate the fluxes of Class I sources in the simulations. 



A 


Aa/Av 


Fx/Fh 


A/L 


Sensitivity limit 


pm 






(/xm/Lo) 




2.2 


0.112 


0.979 


0.685 


16.9 


3.5 


0.063 


0.787 


2.107 


14.5 


4.6 


0.048 


0.714 


9.880 


14.0 


5.8 


0.048 


0.645 


14.65 


13.0 


8 


0.048 


0.654 


27.51 


11.5 


24 


0.005 


1.326 


495.0 


6.8 



